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Abstract 

A self-focusing of a coasting relativistic beam in a plasma channel that is confined by an external 
magnetic field is studied as a means of reconditioning the beam emerging from a beam injector 
[a radio frequency quadrupole (RFQ)] for a linac. A detailed study of the beam stability in the 
self-focused beam has been carried out. In order to explain beam filaments and the resistive hose 
instability in a unified way, we treat all the azimuthal modes in the derivation of the dispersion 
relation in a finite plasma channel that exhibit many unstable modes, which are classified by 
Weinberg's scheme [Steven Weinberg, J. Math. 8, 614 (1967)]. To overcome the energy requirement 
of a beam injector for a high-current, high energy linac, we suggest to add an energy booster of 
a compact synchrotron to the RFQ. The analysis is then applicable to the charged particle beam 
transport in a proton accelerator, such as Large Hadron Collider (LHC) at CERN or APT at 
LANL. 
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I. INTRODUCTION 



The application of electron beam transport in a cold plasma has been slow, but has made 
great strides in recent years for a beam focusing device [Bj, laser-guided beam transport 
problems [?], and heavy ion fusion [8J. In these applications, the m = 1 hose instability, if it 
occurs, displaces the beam position from the center of beam propagation. In both a beam 
focusing device and heavy ion fusion, this can lead to the displacement of the focal point, 
which may cause the loss of the beam to the wall of the device. Furthermore, instabilities 
with higher modes (e.g., m ^> 1) can lead to beam filament. 

Weinberg [2U ES] studied a relativistic beam instability with perturbations on a particle 
orbit which treats all m (azimuthal quantum number) values to account for beam deforma- 
tion, and examines various types of unstable modes in the entire frequency range. 

Resistive instabilities of a charged particle beam penetrating into a plasma channel were 
studied first by Longmire pp. They are complex phenomena and are still not well understood 
for realistic situations. The hose (m = 1 mode) instability was studied by a number of 
authors [UI2JE]- They have shown that a relativistic electron beam (REB) penetrating into 
a cold plasma becomes unstable to the m = 1 mode. However, their analysis remains valid 
for a background plasma whose skin depth (c 2 /A-kuuj) 1 / 2 much larger than the beam radius. 

Since then a significant progress has been made, the REB propagating in neutral plasmas 
has been studied extensively by many authors jl] . Because of the mathematical complexity of 
the problem, it has been necessary for many authors to adopt simple models which are valid 
for specific problems, so that the applicability of models remains limited. In general there 
are serious difficulties in comparing theoretical analysis of an instability with experiments 
(e.g., either a z-pinch or a plasma produced by a gas discharge, which is then confined by 
a magnetic field). We incorporate a feasible experimental condition of interest in the theory 
and compare with experiment. In particular, we are interested in the feasibility of applying 
the theory of a beam-focusing (i.e., reconditioning) in a plasma channel, which is expected 
to be far more effective than the conventional quadrupole magnet focusing in accelerator 
technology. The basic idea of the technique involves the self-pinching of a high-current of 
charged particles by self-magnetic fields [5]. However, the self-pinched beam is subject to 
numerous instabilities associated with particle motion as described below. 

The previous analysis of beam instability was carried out with a model of a beam pene- 
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trating into a cold plasma which is infinite in spatial extent for mathematical convenience. 
However, this approach can not be applied to the self-focusing problems for which the back- 
ground plasma must be finite. This permits one to impose the boundary condition and 
thereby to obtain the dispersion relation that gives qualitative information on the stability 
of the REB. Moreover, it does not apply to the beam focusing by a plasmas in a drift tube 
confined by an external axial magnetic field in accelerator technology, nor does it tell how 
the boundaries of a finite plasma channel are formed. 

A similar situation occurs in an ion-beam transport in the accelerator technology. The 
excessive space-charge is the origin of the beam divergence in a high-current, high-energy 
proton linac. One important new problem arises in this connection, however, is how do 
we overcome the space-charge effect in the absence of a energy booster before injecting 
the beam into a linac, which has not been addressed in the previous work. The principal 
difficulty encountered in the development of the high-current proton accelerator is that a 
high-current beam must reach a critical energy to overcome the space-charge effect. This 
difficulty occurs in a situation for any high-current, high-energy proton accelerator with the 
direct application of a radio frequency quadrupole (RFQ) as a beam injector [TQl [13j fT5] . 

The RFQ by Kapchinski and Teplyakov [13] makes use of a strong focusing with rf- 
electrical field, based on the same principle as is used in a quadrupole mass spectrometer, 
and is basically a homogeneous transport channel with additional acceleration. By the 
geometrical modulation of quadrupole electrodes, one generates the axial field. Thus the 
RFQ has a linac structure which accelerates and focuses the beam with the same rf fields. 
However, the very fact that such axial acceleration is feasible by the geometrical modulation 
brings with it new difficulties which seem to be formidable. 

Yet the acceleration by the RFQ is independent of beam velocity with a large radial 
acceptance, which is a great advantage in the design of a proton linac. Thus the devise 
offers the possibility of utilizing it as a beam injector for a high-current linac. In order to 
make use of it as an alternative to dc acceleration where a beam injector is required, it is 
necessary to demonstrate the possibility of reaching energies to about 7Mev in the RFQ 
before injecting the proton beam for a high-current (80mA) linac such as Large Hadron 
Collider (LHC) at CERN as shown in Figure [l| This prerequisite is determined by solving 
the equations of motion in the presence of self-field in a quadrupole structure in a linac 
[H] . The energy requirements impose additional restrictions on the feasibility of the LHC at 
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CERN. 

The requirements have been appreciated for some time now because of our inability to 
invent a new device that overcome the beam energy requirement for a beam injector for a 
high-current, high-energy accelerator; however it was not difficult to overcome the limitation 
for a low-current accelerator such as the proton accelerator at Los Alamos (LAMPF). 

From a theoretical point of view, the reason for this limitation can be understood, to a 
great extent, with the observations that the beam focusing force by a quadrupole magnet in 
a linac is globally a second-order effect, and the geometrical requirements of a high focusing 
and a large magnetic aperture cannot be achieved simultaneously for a diverging high-current 
proton beam (9J [21] . 

There are several review papers on the progress of RFQ [15] and its wide application 
in many laboratories. The maximum attainable beam energy still remains about the same 
in spite of intense effort to overcome the limitation posed by a high-current, high-energy 
accelerator. 

For example, we give here the numbers for a figure of merits for the LHC planned at CERN 
[19] ; with 100mA input to the RFQ, 80mA proton beam was accelerated to 520KeV by RFQ, 
of which 65mA proton beam was accelerated through the first linac tank. Subsequent linac 
tank cannot accelerate the high current beam and loses a significant amount of the beam 
as shown in Fig. 1 which was obtained by solving Eqs. 13 a-b of Courant, Livingston and 
Snyder [§]. It is this fundamental limitation of RFQ that makes it exceedingly difficult to 
build a high-current, high-energy accelerator such as LHC at CERN. In general, the higher 
the current density of a proton beam, the higher the injection energy will be in order for 
quadrupole magnets in a proton linac to transport the beam with a sufficient focusing force. 

The calculation of the above energy requirement is most simply carried out by introducing 
the condition that the particles remain in oscillating orbits in one direction in the presence 
of self- fields and quadrupole magnetic fields with a given field gradient: 

+ (K 2 - Kl)y = 0, (1) 

where (K 2 - Kf) > 0, K 2 = (dB z /dy)/(BR), K\ = J^ff with 7 = 1.0, dBjdy = 
3.94 x 10 3 Gauss/cm [9J, and P ± = 3.0 x 10~* BR(Gauss -cm) [22] • 

The Figure [T] shows the critical beam injection energy for a given current density, and 
shows the domain of beam stability. The occurrence of beam divergence due to the space 
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charge effect is the origin of the beam divergence which results from our attempt to accelerate 
a high-current beam in a linac, since the focusing force by a quadrupole magnet is not 
sufficiently strong to overcome the space-charge, which poses the fundamental limitation of 
the current accelerator technology [Til IT5] . 

Although the demonstration of successful operation of RFQ in low-current proton accel- 
erator at Los Alamos (LAMPF) has raised a hope of developing a high-energy, high-current 
proton accelerator, the fundamental difficulty in the development of a high-energy, high- 
current accelerator is remains unsolved; that the high-current beam diverges (bursts) as it 
emerges from RFQ in a linac is a real challenge to overcome. This difficulty occurs even if 
the best possible design parameters of RFQ to overcome the Coulomb repulsion in the beam 
are employed [15] . Yet the problem of space-charge of a high-current beam was the original 
motivation of Kapchinski and Teplyakov's concept of RFQ as a possible beam injector for a 
linac. 

One important new possibility that arises in this connection, however, is a utilization of 
self-focusing by a neutral plasma channel which has not been treated in the previous papers. 
The existence of the current limitation in the RFQ has made it necessary to study the 
possibility of space-charge compensated beam transport [16J in the development of a high- 
current, high-energy proton accelerator such as a large hadron collider (LHC) at CERN. 

An important question one should ask: is there any method that can be applied to 
reconditioning the perturbed beam emerging from RFQ before injecting it to a linac? A 
well-known application of space-charge compensated beam [TB] is an effective, useful concept 
to overcome the problem of space-charge effect which is the major stumbling block in a high- 
current beam transport. The decisive advantage of this approach was demonstrated in a 
high-current electron accelerator by shielding the space-charge with a quiescent plasma [6]. 
The focusing force may reach to the value, that is greater than that of a super-conducting 
magnet by several order. Here we assume the presence of an over-dense plasma in which the 
plasma density is much higher than that of a beam. Thus the Coulomb repulsive force due 
to space-charge in the beam is balanced by the self-field of a bunched beam maintaining a 
constant radius. However, the effectiveness of the space-charge neutralized beam transport 
depends on (3 = v/c. Consequently, since the beam emerging from RFQ is in non-relativistic 
domain ((3 = v/c <C), the feasibility of a beam reconditioning is out of the question. This 
technique can be applied only to a charged particle beam with the velocity f3 — v/c > 1 by 
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the combined use of the synchrotron as the energy booster. 

For some time it has been realized that it might not be possible to make use of the RFQ 
as a beam injector for a high-current high-energy linac [TJ)J US]- The simplest and probably 
the best, way of accelerating the high-current proton beam is to apply the phase-locking 
method in the synchro-cyclotron (TTJ [18], but it will be difficult to efficiently extract the 
high-density beam. Yet it seemed feasible to attempt the experiment with a combination 
of a synchrotron with the RFQ , instead of the betatron injection [20], to accelerate the 
beam to the relativistic domain, in spite of the fact that the beam is being lost in the 
extraction process. Thus the role of the RFQ is promising as a pre-injector for a high- 
current, high-energy accelerator provided that the beam emerging from the RFQ can be 
reconditioned by a plasma channel described below. In PHERMEX facility at Los Alamos, 
the high-energy 7-rays are produced by short pulsed high-energy electron beam. To make 
a uniformly diverging x-ray source, it is essential to have a highly focused beam, a pencil 
beam, that can be achieved only by a plasma focusing. 

Finally, the dependence of (3 = v/c on the effectiveness of space-charge compensated 
beam transport arises from the form of self magnetic filed which is the primary focusing 
force. Since the stability due to beam modulation has been studied in a plasma channel 
[21] and a somewhat more detailed study of the model has been already made by Weinberg 
himself [25J, here we study the beam reconditioning with the assumption that the beam 
is accelerated somehow to the relativistic domain in which there is no beam modulation 
and is coasting in a plasma channel confined by an external magnetic field. In this process 
the beam is being reconditioned and is to be transferred to a linac in the next stage of 
acceleration [19J. As we shall see, this simple model calculation presents a rather formidable 
mathematical challenge which defines the ultimate fate of LHC at CERN. 

We shall be interested here primarily in the stability of the relativistic beam, in particular, 
the hose instability of a coasting relativistic beam in a plasma channel confined by the 
external magnetic field since this is the one which may lead to the loss of charged particles 
before the beam is injected to the main linac. An elegant formalism of beam instability is 
provided by the use of the Lagrangian displacement vectors which have been introduced by 
Bernstein, Frieman, Kruskal and Kulsrud in their study of a Rayleigh type energy principle 
for hydro-magnetic stability problems [29J and is applied to the Boltzman-Vlasov equation 
for plasmas by Low [32] . The treatment given by these authors has shown that the method 
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is indeed very powerful in treating stability problems. Here we make use of the concept of 
Lagrangian displacement vectors |29j to give an alternative, more mathematically elegant 
method leading to Weinberg's results |25J in the limit of vanishing external magnetic field 
[36]. 

The present model is essentially equivalent to that of Weinberg (25] , but is different in one 
major respect, that is, we incorporate the external magnetic field to confine the necessary 
plasma inside the drift tube, so that we may impose proper boundary conditions. The 
space-charge compensated beam transport is essentially the same for both the high-current 
electron rf linac (ILC) and the high-current proton linac (LHC), although it is far more 
useful to overcome the difficulty of focusing the beam in a proton linac. 

With these qualitative remarks as an introduction, we proceed to the development of the 
space-charge compensated beam transport in linac. We therefore focus on the stability of 
the REB penetrating into a plasma channel confined by a weak external axial magnetic field 
(i.e., Bq x1 = 1.0KG), consistent with the experiment [BJ- Since we introduce perturbations 
in particle orbits which are determined by the external magnetic field and the self-field, it 
is important to include all the relevant fields in the equations of motion of a particle. 

It is the purpose of this paper to present the essence of our results and some of the details 
of how the dispersion relation is obtained. The present work is a generalization of Weinberg's 



work on the general resistive instability [25] . To make the discussion of our model clearer, we 




derive the dispersion relation by Weinberg [23] by taking the limit uo c = which implies that 




no external magnetic field applied to the system. This will illustrate the versatility of the 
formal application of the Lagrangian displacement vectors in the study of beam instability. 
It is shown here that, even in the presence of an axial guide magnetic field, the REB remains 
unstable to various perturbations. These include the m = 1 hose instability and the set of 
the higher modes that are unstable to various perturbations (type A, B, A, B, C, D, and 



II. BASIC EQUATIONS 

There are two avenues of approach to this problem. One is to pursue the approach of 
Weinberg [25] and consider the the first-order effects by the external magnetic field. Though 
some details of his theory are open to criticism, there is no doubt that the paper points to 
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the right direction for this difficult problem. However I must confess that I have been unable 
to follow the early part, "Introduction", of his paper. A brief study of Weinberg's paper has 
convinced me to take a more formal approach that is transparent. Indeed it is much more 
preferable mathematically to apply a straightforward technique to a complex perturbation 
problem of a many-body system. Moreover, without a confining external magnetic field, 
Weinberg's model does not allow one to impose the necessary boundary conditions that are 
essential to derive his dispersion relations, Eq. 1-9 [25J and is therefore inconsistent with 
his object of the problem that he set out to solve. Hence the object of this paper is also to 
establish what is correct in Weinberg's analysis, and to remove some of ambiguities. 

To obtain a proper dispersion relation, it is first necessary to specify the unperturbed 
motion of an electron in the beam. The REB is assumed to move in the z-direction with 
the average velocity v. The self- field is then given by 



where p is the hydrodynamic pressure on the beam. The solution of Eq. d3j) is then 




(2) 



where n(r) is the particle density of the REB. 

Then the pinch effect can be described by the force (MHD) equation 



(3) 



dr dr \8tt J 47rr' 




(4) 



with the axial pressure po 




(5) 



Hence Eq. Q can be written 




(6) 



In addition we define 




(7) 



which is equal to a 2 

4ne 2 n/('jM). 



e(3/('jM)(B s /r) : where (r 2 a 2 )' 
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>p and ujp 
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The rotational motion of a coasting beam can be described by the equation of motion in 
the presence of an external and a self magnetic fields: 

(- + v ■ V)(v e e e + v z e z ) = -^(/3 x B ext + (3 x B s ), (8) 

where uo c = eB ext /(jMc) = ecB ext /S and v z = constant. 

One can easily show from Eq. d8| that the angular frequencies are given by 



u;± = ^[-o; c ± V^ c 2 + 4a 2 ]. (9) 

Hence the transverse motion of a particle in the polar coordinates can be described by 
the equation: 

r = r x uj c — a 2 r, (10) 

where u c = eB ext /( / yMc). Here 7 and M are the Lorentz factor and the mass of a charged 
particle respectively. 

If we now define the rotational velocities as u± = r — u> T e z x r, then we can show that 
(uj + — uj-)r = —e z x (u + — u_), which describes a uniform circular motion of an electron with 
frequencies u± in the r — 6 plane. Furthermore, by taking the Landau gauge A = ^B^ xt x r, 
the canonical momentum becomes p = / yM(u + + u_)/2. The rotational velocities are not 
a constant of motion; but the canonical momentum p can be related to the Hamiltonian of 
the system [23] and hence we employ p in the derivation of a dispersion relation. 

If a beam is in a relativistic motion in the presence of electric and magnetic fields, then 
its translational motion of a particle in a fluid element must be in a manifestly covariant 
form [22]: 

f = ^*+i**B-wm, (id 

where f3 = v/c. 

Next for a uniform beam with density profile, 




n(r) = { (12) 

the rotation frequency is given by 

a\r)=l " ' (13) 

I JpC? jr 2 if r > a. 

We often simply write for a uniform beam a 2 = ul = 2ne 2 n(5 2 /7M if there is no ambiguity. 



III. FIRST-ORDER EQUATIONS 



The treatment of the first-order equations of motion in Weinberg's paper did not appear 
to be in the most convenient form which is difficult to follow through [22]. 111 fact if is 
almost impossible to introduce an external magnetic field in a formal analysis in Weinberg's 
approach, which is necessary for a systematic study of the beam instabilities in a finite 
plasma channel. 



A. Calculation of Perturbed Fields 



In this section, first we write down the first-order field equations from the Maxwell equa- 
tions assuming that all variables vary as fi(r) exp(im8 + ikz — iujt): 



q 2 E lr 



q 2 E w 



q 2 B ie 



q 2 B lr 



mu .,dE l2 

B lz + ik^— 

r c or 

m .ujOB 1z 

— kE lz - % — 

r c or 



UJ 



A-Ki — J lr 



UJ 

Am — J ig , 
& 



.udE lz m k 

i k—B lz - A7n-J lr , 

c or r c 

muj ., dB lz ,k 

E lz + lk— h 4:711- Jig, 

r c or c 



(14a) 
(14b) 
(14c) 
(14d) 



where q 2 = uj 2 /c 2 — k 2 . 
com 
dE l 



Second from Eqs. ( 14 ) we obtain the decoupled field equations: 
1 d 



r dr 



Or 



m 



E lz + q z E u 



Attuj 



c 2 (k 2 + q 2 ) 



m 



1 d ( dB l2 

I r 

r dr \ dr 



m 



B lz + q 2 B lz 



47T 

rc 



-iq J lz + i — kJw + 







k d 
r dr 



rJir 



(15a) 



(15b) 



Here we shall examine in detail the first-order equations of motion by introducing the 
displacement vector in a particle orbit £(r , t), where Tq describes the unperturbed trajectory 
of a charged particle (an electron or a proton) and t is the time. Here the displacement £ is 
defined by the equation r = r + £(t"0)^)- Upon introducing the Lagrangian displacement 
vector it is possible to expand the velocity in terms of This makes it simple to derive 
the first-order equations of motion in the presence of an external magnetic field and the 
perturbed fields. In particular, we must modify the equation of motion for a fluid element 
that moves with relativistic speed in the presence of electric and magnetic fields as described 
in Eq. (fill. 
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Next we expand the velocity to the first-order in £ defined by r = vq + limiting to the 
fast mode for the time being, 



v{r + £) = «o + «i(£)> 



(16) 



where vi is given by the following equation. 



(17) 



where we have limited to the fast mode u + in Eq. ^ and will repeat the same calculations 
later for the slow mode u- to check algebras. 



We then derive the first-order equation of motion from Eqs. (11) and (14) after a short 
algebra: 



_| + 2vo ■ + (« • V)(« • V)£ = -[£i + /3x Bi + ftx 

+/3x x B S + /3 x (£ • V)B S - (3((3 ■ Eh.)], 



(18) 



where (3 = v/c and the subscript 1 denotes the first-order. 



B. Calculation of Perturbed Density 

To calculate the perturbed current density, we must still find a way of expressing the 
perturbed beam density in terms of displacement vectors. A simple alternative method to 
that of Weinberg [25] is to linearize the equation of continuity: 

— + V(n«) = 0, (19) 

with care on x — > x + £ and n(x + £, t) = no + (£ • V)no + rii and then to pick up the 



first-order terms in the expansion of Eq. (19): 

_1 + V • (n «i) + V ■ (th« ) + («o ■ V)(£ ■ V)n - V(£ ■ V) • (no« ) = 0, (20) 

where we have used V — >• Vo — Vo $, ■ Vo in the expansion and then dropped the subscript 
from Vq. 



Returning to Eq. (20) with ^ = ^ i ^ e l ( kz + md anc l expressing £ r and we obtain 

— — mo; + )ni — i(fi — muj + ) V ■ (n £) = 0, (21) 
11 



which yields the perturbed density n\ = — V • (no£). This is exactly the same as Weinberg's 
derivation of the perturbed density which is much more elegant a method of deriving the 
perturbed density (Eq. (5.10) of Weinberg's [2"5]). 
Hence we write the density perturbation as 

n(x + $) =no(l- V-£), (22) 

for a uniform beam. 

This is perhaps the most crucial equation that describes the collective effects. We can 
impose an appropriate jump condition at the beam-plasma boundary with the equation 
and derive the first-order equations of motion in terms of the displacement vector £ in the 
many-particle system. As mentioned earlier, there is a domain of parameters in which we 
can find a tractable solution for the hose instability in a plasma channel, for which we refer 
the reader to Weinberg' paper [25J. 



IV. DISPERSION RELATIONS FOR A UNIFORM BEAM 



Since the hose instability in a modulated beam that includes the beam bunching in the 
non-relativistic domain is exceedingly complicated, we limit the hose instability in the rela- 
tivistic domain in which a model illustrates a charge-compensated beam transport consistent 
with the experiment [6]. Moreover, we set aside the question of whether a reconditioning of 
a relativistic beam emerging from a beam injector can be studied quantitatively with the 
aid of the coasting beam model in a plasma channel. Hence we leave the question of an 
appropriate beam injector open for future discussion. 

To obtain the desired dispersion relation for a coasting beam for which wake fields are 
negligible, we assume the first-order quantities vary as f ie - t (^ jt ~ m9 ~ kz ) as before. Here we 



limit our calculation to the fast mode uo + . Then Eq. (18) yields after a brief algebra, 

- [(n - rnw+f - r{a 2 )%. + z(fi - muo + ){2uo + + u e )£ g = ^(E lr - (3B W ), (23a) 

- mu + ){2u + + uj c )i r -(n- mu + ) 2 ie = -^t(Ew + 0B lr ), (23b) 

71M 

-(n-muj + )^ z = -^E lz} (23c) 

where Q = to — kv. 
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The significance of the right hand-side of Eq. (23c) is easily seen if we compare it with 



those of Eq. (23a) and Eq. (23b). In the relativistic domain for which 7 ^> 1, Eq. (18) is 
essentially reduced to a two-dimensional (r and 9) problem, since £ 2 is too small to retain and 
may be dropped in the analysis. Hence there is no beam bunching in the relativistic domain. 
This observation is consistent with the order-analysis of Weinberg [241 125] and gives a hint 
that there is a domain in which we may find a tractable solution for the beam instability. 
The straightforward solution of these equations is still very difficult although perhaps not 
impossible. Hence we follow Weinberg's analysis very closely to obtain approximate solutions 
in the limiting cases. 

With the same set of basic assumptions as in Weinberg's analysis \u\a <C 1, \k\a <C 1, 
and \q\a ~ 1, we look for the parameters of the beam and plasma c 2 / (Airavo) <C 1 and 
uipa/vo <C 1 such that ua <C vq and \k\a <C 1. Repeating the same order analysis of J\ r , J\g 



and Ji z in Eq. (15a), we arrive at the equation: 



Id f dE lz \ m 2 2 — Anie(3uj 



r- 



- —E lz + q 2 E lz = y —n x (24) 



In the following analysis we follow closely the procedure of Weinberg [25] and make the 
same basic assumptions, \u\a <C v , \k\a <C 1, and \q\a ~ 1, where q 2 = —k 2 + (iuj/c 2 )(A7Ta — 
ioS) —* Airiau/c 2 . Here a is the beam radius. 

By trivial extension of the arguments leading to Eq. (3.30) of Weinberg [25], we obtain 

--rr—E -—£ + q 2 S = n x , 25 

r dr dr r z cr 

where E lz = g(j^ e - wt + lkz + ime . Here we have taken the plasma current J p = cxEi, where 
cr is a scalar. In a dense plasma, Ohm's law in its simple form remains valid in a wide 
frequency range. Hence the Hall effect in the REB is negligible. This assumption of scalar 
conductivity in our model is reasonable, because io c jv e i 1, where uo c ~ 5.0 x 10 8 sec _1 and 
v ei ~ 2.2 x 10 9 sec -1 . Here a ~ 1.4 x lO^sec" 1 for n p ~ 10 12 /cm 3 at kT = 12eV [22]. 



We rewrite Eqs. (23) to calculate n\ 



- [(0 - muj + f - r(a 2 )% + i{Q - mu+)(2u + + u c )£ 9 = (26a) 



+i(n-muj + )(2uj + + uj c )£ r + (n-muj + ) 2 £ e = (26b) 

7M 



where Vt = uj — kv and (or)' = in Eq. (26a) for a uniform beam 



13 



Here we have extended Weinberg's analysis on the right-hand sides of Eq. (23a) and 



Eq. (23b) to write 



w 



(E\ r — (3B ie ) 

:{E 1Q + j3B lr ) 



tCV J~ 17") 

-E' lz = exp [— i{ut — mO — ikz)\ - 



7Mu 



7 M 



-E\ z = exp [— i(ut — m6 — ikz)]^'^ 



7 AP ' " 7 M /J 7 M ' 

where we take approximate values of Fi in terms of the electric field E\ z [25J. That is 

— iev 



(27a) 
(27b) 



Mr) 



u 

evmS(r) 
uj r 



Solving for £ in Eq. (26|) in terms of the perturbed field E\ z < 
2iev 1 
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S'(r) - {2UJ+ + Uc) (-) £(r) 

jMuE(u)[ K > (fi-mcu+) \ r J K [ 

ev I r (2u; + + u>c 



(28a) 
(28b) 

-iwt-ime+ikzg^ we have 
(29a) 



r{a 2 )' 



(fi — muj + ) 



S'(r) + [1 



(fi — muj + ) 2 \ r 



(29b) 



where H(a>) = (fi — mco> + ) 2 — 4a 2 — r(ct)' — uj 2 



Making use of the identities n(r ) = (r 2 a 2 )'(7M)/(47re 2 /3 2 r) and 2lu + + lu c = y/w* + 



a* 



and defining m = e < m9 + kz w )A/"(r) 
of £ as 



-V • (n£), we obtain the perturbed density in terms 



Af(r) 



ic 



where 



/+ 
9+ 



1 d ( £ \ nr 
r dr V dr 



2a 2 + r(a 2 )' 



f+S - g+£ 



uj 2 + 4a 2 - (fi - muj + ) 2 + r(a 2 )' 
m\ 1 d 



V^+4^/4 



(30) 

(31a) 
(31b) 



r / fi — mu + dr 

For the angular frequency cj_ = |[— c<j c — ^Juj 2 + 4a 2 ], an exactly parallel calculation 
yields /_ and g_ . This result shows an interesting symmetry that /+ — > /_ by means of the 
substitution cj + — ► cj_ or vice versa. It should stressed that this calculation for /_ and g_ 
should be carried out to make sure our algebras are indeed correct, although it is somewhat 
tedious. 

Since the electron motion in the r— 9 plane has slow and fast rotations, and the generalized 
momentum (canonical momentum) is the average of the mechanical momenta in the presence 
of a static magnetic field, we take /=[/+ + /_]/2 which is given by 
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/ = [2a 2 + r(a 



2 VI 



e 2 



(32) 



6 4 - m 2 {uj 2 c + 4a 2 ) (Q + mu c /2) 2 ' 
where 9 2 = (1 - m 2 /A)(uj 2 c + 4a 2 ) - (Q + mu c /2) 2 + r(a 2 )'. 

It should be stressed that an argument similar to that of Weinberg [25] for taking an 
average based on the probability of rotation in either (positive or negative) direction may 
not hold in the presence of an external field. Since the equation of motion can be written as 
dp/dt = e[E + /3 x B], it is indeed correct to take the canonical momentum in our analysis. 
This clearly shows the inadequacy of Weinberg's model. And yet in the limit of B^ xt — > 0, 
our final results reduce to those of Weinberg [25] . 

Similarly, we take g = [g + + g~}/2 and then is given by 



9 



( m 2 
\ r 



(Q + mu c /2) 2 
(Q 2 — m 2 a 2 + muj c Vt) 



(u 2 c +4a 2 )f(r) 

02 



+ - 



m 
r 



2 V r ) Vt 2 — m 2 a 2 + mui c Vt 
It should be noticed that, for a uniform beam, (a 2 )' 



(33) 



which simplifies algebras im- 
mensely. Henceforth we will assume that a uniform beam is coasting in the over-dense 
plasmas and that the beam is completely neutralized. For the special case of a uniform 
beam density profile, it is possible to obtain a wave equation similar to that of Weinberg 
[251. 



By taking n(r) = n6(r — a), it follows from Eqs. (32) and (33) 



f(r) 



1 — rf if r < a, 
if r > a. 



Here rj is defined as 



7] 2 = 1 



Similarly, 



2[(4 - m 2 ) + (1 - m 2 /A)\ 2 c - A 2 ] 
[(4 - m 2 ) + (1 - m 2 /4)A 2 - A 2 ] 2 - m 2 (4 + A 2 ) A 2 

9(r) = -(- 



where 



C 2 = 2m 2 (l -rf 



1 + A 2 /4 



A 2 - m 2 (l + A 2 /4) 



1 + 



2A 2 



A 2 -(l-m 2 /4)A 2 -(4-m 2 ) 



(34) 



(35) 



(36) 



(37) 
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with A = A + mA c /2, A 
uniform beam density. 



fl/a, and A c = u> c /a. Here a 1 



27ce 2 nf3 2 /{-fM) for a 



Substitution of these expressions for the perturbed density in Eq. (25) yields 



1 d 

r dr 



.dS\ 



m 



r[l - 1 ~ ^[1 - f(r)}£ + [q 2 + g(r)}£ = 0. 



(38) 



The detailed algebra leading to this wave equation is straightforward, but it is somewhat 
tedious. Similar to that of a two-dimensional vibrating membrane problem, the solution of 



Eq. (38), which is finite and satisfies the boundary conditions for a finite plasma channel 



is given by Hankel functions 



£(r) 



Here H^(qr) = J m (qr) + iN m {qr) 



J m {qr/rj) if r < a, 

H L(l r ) ~ a mJ m {qr) if r > a. 



(39) 



8{r) by Eq. (39), together with the boundary conditions at the surface of the plasma 



channel of radius b specifies the eigenvalue a m , 

#i-i(?&)/^m-i(<Z&) ifm>l 
H^\qb) I 'Ji(gfe) ifm = 0. 



a:, 



(40) 



Integrating Eq. (38) over a — e and a + e, we obtain 



a£\a + e) - arf£\a - e) = £ £(a) 



(41) 



Substitution of appropriate solutions of the wave equation from Eq. (|39|) into the left hand 

or 

c 2 



side of Eq. (41 ) and rearrangement of terms in the limit e — > yields the dispersion relation: 

[H^\qa)]'-a m {qb)J' m {qa) 



J m (qa/r))J qa 



Hln\qa) - a m (qb)J m (qa) 



(42) 



where m > 0. 

This dispersion relation is identical in form to that of Weinberg [25] , but differs in r\ and 
(. In the limit of vanishing B ext , the dispersion relation goes over into that of Weinberg [25J. 
But it should be noticed that the boundary condition at the edge of the plasma column in 
Weinberg's analysis [25] is not valid, since there is no external magnetic field that confines 
the plasma channel in experiments. It is therefore apparent that his dispersion relation 
Eq. (1.9) [and Eq. (12)] is inconsistent with the problem he has posed in his paper and is in 
self-contradiction. 
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V. ANALYSIS OF THE DISPERSION RELATION FOR A UNIFORM BEAM 



We shall be interested here mainly in the resistive hose (m = 1) instability in the low and 
high frequency limits since the resistive hose mode is the one that leads to the loss of a beam. 



We note from Eq. (42) that, since q = Aniau/c , it may be possible to obtain a rather simple 
solution which holds to a higher order of approximation in the low and high frequency limits. 
The method of making such an approximation lies in the realization that the conductivity 
of the plasma remains fixed for a given plasma density and the approximation of Bessel 
functions in the asymptotic limits as a function of us readily available. 



Thus the analytic solutions of Eq. (42) can be studied in low and high frequency regimes 
with the asymptotic limits of Bessel functions; the two asymptotic solutions should then be 
connected smoothly by analytic continuation. For each m, we may classify the modes as 
in Weinberg [25]. This classification is not entirely trivial, since the external magnetic field 
introduces new modes by removing the degeneracy found in Weinberg's analysis [25]. The 
central question in determining the efficiency of beam transport by means of plasma focusing 
is to find which instability would affect the beam transport most significantly. We present 
here only an outline of the classification similar to that of Weinberg [25] with emphasis on 
the resistive hose mode (|u;| <C a) that affects the beam transport most dangerously, since 
if the m = 1 hose mode occurs, the entire beam can be lost. Moreover the two-stream mode 
for which uo ~ o has been already treated in detail 



Returning to Eq. (42) we expand the left-hand side (LHS) and the right-hand side (RHS), 



using the following identities in Bessel functions: 



H${qa) = Jm{qa) + iY m (qa), 

J' m {qa) = (l/g)[gJ m _i(ga) - (m/a)J m (qa)), 
ffffV) = (l/q^qHUUqa) - (m/a)H^(qa)], 



as 



LHS = Tj 



J m -i(qa/ri) r) , 

— ffiy — 

Jm{qa/v) Q a ' 



(44a) 



RHS _™ + Vm-iCgo) J m _i(g&) - y m _i(g6) J m _i(ga) 
qa Y m (qa)J m ^i(qb) - K m _i(g6)J m (ga) 
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With Eqs. ((44aJ) and ( |44b[ ), Eq. (|42]) becomes 
J m -i(qa/r)) 



T] 



J m (qa/if) 



m 



--(i-v 2 )-- 

qa qa 

Y m -i(qa)J m -i(qb) - F m _i(gb) J m _i(gg) 
Y^(ga) J m -i(qb) - Y m _i(qb)J m {qa) 



(45) 



Here if and C 2 are defined by Eq. (35) and Eq. (37). 



A. Low Frequency Regime: \q\ — ► 



For | go | <C 1, the dispersion relation Eq. (42) can be rewritten as 



V_ ( J' m {(i a h) \ 

qa \ J m (qa/rf) J 



(qa) 2 + 1 h 



(46) 



where m^O, and 



(47) 



where m = 0, qb <C 1. 

Here we have used the identity Hm(qa) = J m (qa) +iY m (qa). We now classify the modes 
into those for which \qa/r)\ — > at low frequency and designate C and .D modes, otherwise 
we call A and I? modes. 



For m = 1, ga C 1, and g& <C 1 in Eq. (45), we obtain the following expansion of the left 
hand-side of the equation: 

r\ J[(qa/r)) if (qa) 2 <? 4fl4 

qaJi(qa/rf) (qa) 2 (^v) 2 9Gf] 4 

and similarly the right hand-side can be expanded as 

1_ Y (qa)J (qb) - Y (qb)J (qa) 

qa Y 1 (qa)J (qb) - Y (qb)J 1 (qa) 



1 

ga 



1 - (qa) 2 \n(b/a) + ^(ga) 4 ln 2 (6/a) + iln(6/o) 



1 

ga 



(ga) 4 (6 2 - a 2 ) + 



Next we rewrite Eq. (35) as 
rf = 1 - 



[3(1 + A c 2 /4) - A 2 



[A 2 - 3(1 + A 2 /4)] 2 - 4(1 + A 2 /4)A 2 

l-A 2 /2 + ---+)A 2 + ■■■+, 



1 , ^.M 
3 6 c 27 



(48) 
(49) 
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and similarly rewrite Eq. (37) 

c 2 



2(1 + A 2 / '4) 



*»-(! + XI/ A) 



1 - 



2A 2 



3(1 + A 2 /4) - X 2 



(50) 
(51) 



where a 2 = A 2 /(l + A 2 /4). 



As shown by Weinberg [2S], for \qa/r]\ C 1 we substitute the expressions Eq. (48) and 



Eq. ( 50 ) to obtain an implicit dispersion relations 

1 



A + A H h 



4 A 2 (1 - A 2 /2) - 1(1 - A 2 /2){i(l - A 2 /4) + ln(6/a)}(ga) s 
+ \{l ~ Ae){^ + ^A 2 (l - A 2 ) + ln 2 (6/a) + ln(6/a)}( ? a) 4 
- hi- XDqW + ■■■ + . 



(52) 



Here we have assumed \qa\ -C 1 and \qb\ <C 1. In the limit q 2 — > 0, A is oscillatory, which 
is the hose instability and shows that, in the presence of the external field, the growth rate 
is reduced due to the restoring force by the magnetic field. This is in agreement with the 
numerical result although the effects are not significant [37]. 

While it is possible to calculate the A 2 to the order of (ga) 4 by an iteration technique as 
in Weinberg [25J, we just write, for the sake of simplicity, A 2 to the order of q 2 a 2 instead: 



A + A C A 



-A 



A 2 /2){^(l-A 2 /4) + ln(Va)}M 2 + 



+ 



(53) 



Here A = A + A c /2, A = ±(u — kv)/up, and q 2 = (Aniau) / c 2 , and A 2 = uj 2 /uj 2 . 
Thus, 



A 



-Ac ± (1 - A 2 ) 1 / 2 {A 2 - 2[\(1 - A 2 /4) + ln(6/a)]g 2 a 2 } 1/2 



(54) 



Eq. (54) is the dispersion relation for the hose instability in the limit q a — > and it 



shows that, in the presence of the external magnetic field, the growth rate is reduced due to 
the restoring force by the magnetic field which is in agreement with numerical computations 
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B. Unstable Modes (A and B): \qa/rj\ does not converge to zero: 



If qa/rj does not converge to zero, rj must then converge to zero as fast as \qa\. Rewriting 



Eq. (35) as 



= A 2 -2A-m 2 (4 + Ag)A 2 
' A 2 -m 2 (4 + A 2 )A 2 ' 1 ' 

where A = (4 - m 2 ) + (1 - m 2 /4)A 2 - A 2 . 

Hence the numerator of the above equation must be zero which defines the modes of A 
and B modes: 

_ o f [(1 + m 2 /4)A 2 + (3 + m 2 )] - A^ 2 (A mode), 

A = < (56) 

[ [(1 + m 2 /4)A 2 + (3 + m 2 )] + A\L 2 (B mode), 

where A m = 1 + 4m 2 (l + A 2 /4)(3 + A 2 ). 

With the external magnetic field, the degeneracy of m = and m = 2 are removed in A 
and B modes found in Weinberg's analysis [25] . 



But as in Weinberg's analysis, m + ( 2 in the righthand side of Eq. (46) does not vanish 
for m > when A 2 takes either value of the above A and B modes. 

Hence we may carry out similar analysis iteratively for A 2 from Eq. (46) for A mn - mode 
by the equation 

V -» (qa)/jmn, (57) 
where j mn is the n-th root of J m (x) = which yields 

:a ,o , . , _2^a .fW 1 M 2 ro , m 2 (4 + A 2 ) - 2 1 



Xi^ (3 + m 2 ) + (1 + m 2 )A 2 -VA 1 --^[2+ v ' cJ ]. (58) 

Z Jmn V ^1 



Similarly, for 5 mn - mode we have 

p _ (3 + m >) + (1 + m 2 )A? + VA 1 + i%l![2 + m2(4+ ^»- 2 ], (59) 

2 Jmn V Ai 

where Ai = m 2 (4 + A 2 ) (3 + A 2 ) + 1. 



Next for m = 0, Eq. (47) can be written as 

rj^(qa)/y n . (60) 
Here y n is the n-th root of the equation 

J'MKvHy)) = I {b 2 /* 2 - 1) • (61) 
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Hence for m = 0, there is no difference between the present analysis and that of Weinberg's 



analysis except for the presence A 2 term which removes the degeneracy. Solving Eq. (55) by 
iteration, we obtain 

A 2 = A c 2 + 3 ± 1 - 2{qa) 2 /y 2 n + ■■■ . (62) 



C. m = 1 Resistive Hose Instability: (C, D) modes 



Suppose qa/rj — ► 0. For the left-hand side of Eq. (46) then becomes 

-.2 



v ( JL(w/v)\ _ v 

— — ; j-r ~ TTl — 



qa \J m (qa/r]) J q 2 a 



2 n 1 ' 



(63) 



with the aid of the identity 



J'miia/v) = -J m +i(qa/v) + (mri/qa)J m (qa/r}). 



And yet the condition that qa/r] — > implies that, by Eq. (35), A does not take the 
following value 



A 2 = (1 + m 2 /4)(4 + A 2 ) - 1 ± [m 2 (4 + A 2 )(4 + A 2 ) + 1] 1/2 , 



(64) 



for which rj = 0. 



With the condition \qa/r)\ — >• and with Eq. (63), Eq. (46) yields the following condition: 



mrj 2 + m + £ 2 = 0, 



(65) 



where m > 0. 



The zeros of this equation in which rj 2 and ( 2 are defined by Eq. (48) and Eq. (50) defines 
the modes of type C and D and is given by: 

A 2 = (m 2 -2m+3/2)(l+A 2 /4)±(l+A 2 /4) 1 / 2 [{3m(m - 2) + 9/4}(l + A 2 /4) - m(m - 2)] 1/2 , 

(66) 

where ± signs correspond to the modes of type C and D respectively based on Weinberg's 



classification [25]. Note also that if A — > A for A c = uj c /a = 0, the results of Eq. (66) in the 
limits go over B14 and B15 of Weinberg (25] as they should. 
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D. High Frequency Regime: \q\ — > oo 



1. A Mode and B Mode 



For \qa\ 3> 1 and \qb\ 3> 1, with the aid of asymptotic formula of the Bessel functions of 
the first kind J m _i(ga) and the second kind F m _i(ga) for qa ^> 1, 

o -|l/2 



and 



F m _x(ga) 



J m _i(ga) 



7rga 



cos[ga - ^ - -(m - 1)tt], 



■,1/2 



sin[ga - | - ^(m - 1)tt], 



the dispersion relation Eq. (45) reduces to 



J; 



(67) 



(68) 



(69) 



.J m {qa/v)J 1 a Q a ' 
Next we divide the modes into those for which l^a/r/l remains finite with \q\ — > oo (^4 
and S modes) and those for which \qa/rj\ — > oo {T> mode). 



With some algebraic rearrangement, we may rewrite Eq. (35) as 

2[(l-m 2 /4)(4 + A 2 )-A 2 ] 



V 



[A 2 - (1 + m/2) 2 (4 + A 2 )] [A 2 - (1 - m/2) 2 (4 + A 2 )] 



(70) 



It is easy to see from Eq. (55) that \qa/rj\ may remain finite if the following conditions 
are met: 

1 + m/2) 2 (4 + A 2 ) A mode, 

(71) 

1 - m/2) 2 (4 + A 2 ) B mode. 
For a finite value of \x\ = \qa/rj\, we may rewrite Eq. (69) to obtain the 1/g 2 term in A 2 , 



a! 



JLW + c 2 



xJ m (x) (xr]) 2 



0. 



(72) 



where x = qa/rj. 



The A and B modes for which A^ takes (1 ± m/2) 2 (4 + A 2 ), and using those in Eq. (70) 



and Eq. ( 50 ) we may write 



r±(m, A±) 



(73) 
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where F± is defined as the following: 

2 ' 4(l + A 2 /4) 



r±(m, A±) 



m 



2 VA| -m 2 (l + A 2 /4) 
After a brief algebra, we obtain 



1 + 



A| - (1 - m 2 /4)A 2 - 4(1 - m 2 /4) 



(74) 



r ± (m, A±) = ± 



3 (m ± 2) 
4(m± 1) 



Since I/77 2 — > as — > 00, 



£ 3(m±2) 

n 2 + 



77- 4 (mil)' 

which is independent of A c . If A+ = (1 + m/2) 2 (4 + A 2 ) for A mode, then 

C 2 _ 3 (m + 2) 
r^ 2 4 (m + 1) 

Hence the dispersion relation for A mode is 



(75) 



(76) 



(77) 



A 2 — > 4(1 + m/2) 2 (l + A 2 /4) + 



3(m + 2) x 2 
4(m + 1) (ga) 2 ' 



(78) 



where x is the solution of the equation 



x f (:r) = ± 3 (2 + m) 

± 4 (l + m) ml j ' 

When A 2 = (1 — m/2) 2 (4 + A 2 ) for B mode except for m — 1 and m = 2, we have 



C 2 _ 3 (m-2) 
r? 2 4 (m — 1) 

Hence the dispersion relation for B mode is 

A^4(l-m/2) 2 (l + A 2 /4) + 
where x is the solution of the equation 



(79) 



3(m-2) x 2 
4(m-l) (ga) 2 ' 



(80) 



(81) 



X.J m (x} i t/ m (x), 



(82) 



"4 (1 - m) 

where both m / 1 and m^2. 

We also notice that our results do not go over to those of Weinberg's analysis in the limit 
A c — > 0, because there was an algebraic error in his analysis [see equations (B26) and (B27)]. 
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2. V Mode 



If f) 2 remains finite but ( 2 — > oo, then the mode T> would take place in the beam. This 
happens if A 2 7^ (1 ± /2) 2 (4 ± A 2 ) and A 2 — > m 2 (l + A 2 ) except for m — 1 for which the 
numerator of ( 2 becomes zero with A 2 = m 2 (l + A 2 /4). For m = 0, ( 2 is identically zero. 
Hence the type T> mode begins with m > 2. 

In a high frequency limit, \q\ — > 00, \qa\ ^> 1, and \qb\ ^> 1, T> mode is given by 

p -' 1 + f) + W' (83) 

where x is a root of xJ[(x) = — 3Ji(x) and qa/\r}\ — > 00. 

We have carried out the above analysis to guide numerical work of solving the dispersion 



relation for various instabilities, Eq. (45), by a computer. The above analytical results 
have been borne out in our detail numerical calculations. We have found that the external 
magnetic field reduces the growth rate somewhat, but not significantly since uj c = ecB ext /S 
in the relativistic domain. The results (Figure [2]) are essentially same as the Weinberg's 
analysis. 



VI. DISCUSSION AND CONCLUSION 



In conclusion we note that it may be somewhat confusing by studying both a proton beam 
for a proton linac and an electron beam for PHERMEX facility together, but the stability 
analysis remains valid for both cases with a proper change of charge and mass of a particle. 
The proposed technique of combination of the synchrotron and the RFQ to accelerate a 
high-current beam to meet the necessary energy requirements for the beam injection may 
become feasible, since it takes so short a time to travel for the beam in a short plasma channel 
that the resistive hose instability may not develop to disrupt the integrity of the beam. 
Moreover, the technique has been demonstrated in an electron accelerator by Nakanishi 
group [BJ. However, the experimental evidence we could find (including the plasma cone in 
a Be chamber of PHERMEX facility) is too scant to permit any safe generalization. Indeed 
it would be desirable to have experimental demonstrations for a high-energy, high-current 
proton beam. Yet the case we have discussed so far is conclusive enough to demonstrate 
its feasibility. Without a proper reconditioning a high-current beam by a dense plasma 
channel, a reliable operation of a high-current, high-energy accelerator is highly unlikely. 

24 



However, the creation of a dense plasma-channel may complicate the maintenance of a linac 
in a routine operation and may make it difficult to perform any reliable experiments. Still 
a potential difficulty in developing a reliable beam injector for a high-current, high-energy 
accelerator remains the major stumbling block for the the large hadron collider (LHC) at 
CERN or APT project at LANL. In any practical, realistic sense, the LHC is doomed to 
failure [38j E2J SO] , but particle physicists could still turn to astrophysical observations [41] • 
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FIG. 1: The stable domain, as a function of the beam injection energy and the current, is calculated. 
The condition K 2 - K\ > in Eq. ( 1 ) was imposed on the injection beam energy by K\ = 
for a stable particle orbit. 
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FIG. 2: The growth rate Im(Q/up) of the hose-mode (Ci) as a function of frequency 47ra 2 uj/c 2 
for various external magnetic field strengths The curve for B ex t = 0.0/cG corresponds to the 
Weinberg's solution for which the present model of plasma channel breaks down. As the strength of 
the magnetic field is increased by AB ext = O.bkG, the growth rate decreases, but not significantly. 



29 



